Evolutionary Applications 



Evolutionary Applications ISSN 1752-4571 

ORIGINAL ARTICLE 

Genetic assessment of a summer chum salmon 
metapopulation in recovery 

Maureen p. Small/ Thorn H. Johnson,^ Cherril Bowman^ and Edith Martinez^ 

1 Washington Department of Fish and Wildlife, Molecular Genetics Laboratory, Olympia, WA, USA 

2 Point No Point Treaty Council, Kingston, WA, USA 



Keywords 

conservation, hatchery impacts, 
microsatellites, population structure, recovery, 
salmonids, supplementation 

Correspondence 

Maureen P. Small, Washington Department of 

Fish and Wildlife, Molecular Genetics 

Laboratory, 600 Capitol Way N., Olympia, WA 

98501, USA. 

Tel.: +1 360-902-2682; 

fax: +1 360-902-2943; 

e-mail: smallmps@dfw.wa.gov 

Received: 1 5 October 2012 
Accepted: 16 September 2013 

doi:10.1111/eva. 12118 



Abstract 

Programs to rebuild imperiled wild fish populations often include hatchery-born 
fish derived from wild populations to supplement natural spawner abundance. 
These programs require monitoring to determine their demographic, biological, 
and genetic effects. In 1990s in Washington State, the Summer Chum Salmon 
Conservation Initiative developed a recovery program for the threatened Hood 
Canal summer chum salmon Evolutionarily Significant Unit (ESU) (the meta- 
population) that used in-river spawners (wild fish) for each respective supple- 
mentation broodstock in six tributaries. Returning spawners (wild-born and 
hatchery-born) composed subsequent broodstocks, and tributary-specific supple- 
mentation was limited to three generations. We assessed impacts of the programs 
on neutral genetic diversity in this metapopulation using 16 microsateUite loci 
and a thirty-year dataset spanning before and after supplementation, roughly 
eight generations. Following supplementation, differentiation among subpopula- 
tions decreased (but not significantly) and isolation by distance patterns 
remained unchanged. There was no decline in genetic diversity in wild-born fish, 
but hatchery-born fish sampled in the same spawning areas had significantly 
lower genetic diversity and unequal family representation. Despite potential for 
negative effects from supplementation programs, few were detected in wild-born 
fish. We hypothesize that chum salmon natural history makes them less vulnera- 
ble to negative impacts from hatchery supplementation. 



Introduction 

In the Pacific Northwest of the United States, fisheries 
managers are developing strategies to conserve and 
restore native biodiversity in Pacific salmonids. After dec- 
ades of harvest beyond sustainable rates and loss of 
spawning and rearing habitats, over 20% of native sal- 
mon populations are imperiled or extinct (Augerot and 
Foley 2005). Salmon hatcheries have often been employed 
to mitigate losses of native fish, both as means to support 
native populations and to support harvest opportunities. 
Traditional hatcheries developed broodstocks whose 
progeny were released into nonancestral streams (i.e., 
nonlocal broodstocks); rearing and release practices pro- 
moted adaptation to hatchery environments such that 
nonlocal, hatchery-adapted hatchery fish fared poorly 
under natural conditions and potentially brought mal- 



adaptive traits into wild populations (Tessier at al. 1997; 
Hansen et al. 2000; McClure et al. 2008; Naish et al. 
2008). Currently, impacts from hatchery fish are moni- 
tored by assessing relative fitness, genetic diversity, and 
effective population size (ATe) in supported populations 
(Wang and Ryman 2001; Hare et al. 2011). Recent stud- 
ies have demonstrated lower fitness in populations receiv- 
ing hatchery fish from traditional hatcheries, compared 
with local, wild populations (reviewed in Berejikian and 
Ford 2004; Araki et al. 2008; Christie et al. 2012). Studies 
suggested that hatchery fish with lower fitness or mal- 
adapted traits depressed the productivity of wild stocks 
with which they interbreed (Lynch and O'Hely 2001; 
Araki et al. 2007a, 2008). Further, introgression by com- 
mon farmed broodstock into wild populations may 
decrease genetic distinction among wild populations 
(Hansen et al. 2009; Glover et al. 2012). 
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With concerns over the role of hatchery programs for 
restoring wild populations, supplementation hatchery pro- 
grams were designed to decrease negative hatchery impacts. 
Supplementation hatcheries aim to temporarily increase 
spawner census size to conserve genetic resources and boost 
natural population abundance of imperiled wild popula- 
tions while minimizing genetic and ecological risks com- 
monly associated with traditional hatchery practices (Ford 
2002; Goodman 2004). Supplementation programs devel- 
oped to meet conservation objectives bring local, natural- 
origin adults, juveniles, or eggs into a hatchery to initiate 
the program for a specific tributary (Berejikian et al. 2008; 
Small et al. 2009) and release the hatchery-born juveniles 
into their specific tributary. When adult fish return to the 
supplementation program tributary, hatchery- and wild- 
born spawners intermix in the natural spawning area, and 
subsequent hatchery broodstocks for each target tributary 
are a mix of these hatchery- and wild-born spawners 
returning to their tributary. With supplementation hatch- 
eries, it is hoped that increasing spawner abundance with 
local-origin supplementation fish will have positive effects 
of maintaining genetic diversity, increasing N^, and sup- 
porting persistence and adaptive potential of targeted pop- 
ulations (Wang and Ryman 2001; Hedrick 2005; Berejikian 
et al. 2008). Supplementation programs have met with 
varying success because hatchery practices may counter 
intended benefits through high variance in hatchery family 
sizes and unequal sex ratios, which can decrease and 
reduce genetic diversity of supplemented populations 
(Ryman-Laikre effects, Ryman and Laikre 1991; see Chris- 
tie et al. 2012). Some evidence suggests lower reproductive 
success for natural-origin steelhead (Oncorhynchus mykiss 
Walbaum) whose parents arose from a supplementation- 
style hatchery in comparison with wild-origin steeUiead 
with no hatchery ancestry (Araki et al. 2007b, 2008; Chris- 
tie et al. 2012). Further, Chilcote et al. (2011) suggested 
that these hatchery impacts on wild salmon productivity 
were also found in coho salmon {Oncorhynchus kisutch 
Walbaum) and Chinook salmon {Oncorhynchus tshawyts- 
cha Walbaum) and occurred regardless of number of gen- 
erations in the hatchery or the origin (local or exotic) of 
the hatchery broodstock. However, Sharma et al. (2006) 
observed increased productivity in supplemented wild coho 
salmon and Hess et al. (2012) found higher abundance and 
no fitness loss in supplemented wild Chinook salmon. 

Although chum salmon {Onchorhynchus keta Walbaum) 
have a typical anadromous Pacific salmon life history, neg- 
ative impacts from a supplementation hatchery program 
may be lower for chum salmon than for other salmonids. 
Chum salmon juveniles migrate to estuaries within days of 
emergence (Johnson et al. 1997), and hatchery chum sal- 
mon are typically released as fry after 2-3 months of rear- 
ing (chum salmon juveniles in this supplementation 



program were released into their target streams within 
75 days of hatching). Thus, in contrast to other species, 
hatchery-born chum salmon spend minimal time under 
the selective forces within an artificial rearing environment 
and impacts might be primarily through broodstock selec- 
tion, lack of mate choice, and unequal success of hatchery 
families. In this study, we explore impacts of supplementa- 
tion programs in chum salmon by examining changes in 
genetic diversity and effective population sizes in a 
thirty-year time series starting in the 1970s that includes 
collections spanning the time from before, during, and 
after supplementation. Hatchery supplementation and 
reintroductions began in the 1990s for summer chum sal- 
mon in Hood Canal (HC) and Strait of Juan de Fuca (SJF) 
in Washington State in response to declines and extinctions 
(see Methods for details). 

Materials and methods 

Chum salmon life history 

Chum salmon spawn mainly in lower tributary reaches up 
to barrier falls in coastal tributaries along the northern 
periphery of the Pacific Ocean, as far south as Japan and 
Oregon. Similar to other Pacific salmon, chum salmon are 
semelparous and return to spawn generally after 3 or 
4 years in salt water, but returns after 2 or 5 years are not 
uncommon. In contrast to other Pacific salmon (except 
pink salmon), chum salmon juveniles spend little time in 
freshwater, migrating downstream to estuaries and salt 
water a few weeks after emergence. Chum salmon stray at 
similar rates to other Pacific salmon (Johnson et al. 1997), 
but high spawner densities may increase stray rates 
(reviewed in Quinn 1993). 

Hood Canal summer chum salmon 

Within the Puget Sound region in Washington State, sum- 
mer-run chum salmon in HC and SJF (Fig. 1) are geneti- 
cally and ecologically distinct from fall-run chum salmon 
in the region (Phelps et al. 1994; Tynan 1997) and are des- 
ignated a separate ESU (Johnson et al. 1997). Due to popu- 
lation declines and extinctions, the HC summer chum 
salmon ESU was listed as threatened under the Endangered 
Species Act in 1999 (http://www.nwr.noaa.gov/ESA-Salmon- 
Listings/Salmon-Populations/Chum/Chum-Status-Reviews. 
cfm). A recovery program started in the early 1990s, prior 
to the ESA listing, enacted crucial harvest reductions and 
supplementation programs. In 2000, state and tribal co- 
managers completed the Summer Chum Salmon Conserva- 
tion Initiative (SCSCI; WDFW and PNPTT 2000), a recov- 
ery plan that formalized and expanded recovery efforts. 
Two independent summer chum salmon populations were 
designated as ESU recovery units: HC and SJF; both 
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6. Big Beef Creek if 

7. Seabeck Creek 



8. Stavis Creek 

13. Anderson Creek 

16. Dewatto River 

17. Skokomish River 

18. Tahuya River 

19. IVIission Creek 




Figure 1 Map of Hood Canal, Strait of Juan de Fuca and portions of Puget Sound. Rivers and streams are numbered and asterisks mark tributaries in 
the study. At right are listed the tributary names associated with numbers and the status of subpopulations in the tributaries. Map was modified from 
Sands et al. (2009). 



consisted of multiple subpopulations in rivers draining to 
those two marine basins and both were monitored for sta- 
tus relative to recovery goals. There had been small releases 
of hatchery summer chum salmon juveniles in various HC 
tributaries before 1938 (Johnson et al. 1997), but the recov- 
ery program was the first directed hatchery intervention in 
most tributaries. 

As part of efforts to restore wild subpopulations of sum- 
mer chum salmon, supplementation hatchery programs 
were initiated for three subpopulations (Lilliwaup, QuQ- 
cene, and Salmon) in 1992 and for three more subpopula- 
tions (Hamma Hamma, Jimmycomelately, and Union) in 
1997, 1999, and 2000, respectively (each supplementation 
program used wild-origin broodstock collected from its 
respective tributary). Summer chum salmon from three 
streams with supplementation programs (Quilcene, 



Salmon, and Union) were also the donor stocks for pro- 
grams used to reintroduce summer chum salmon in one 
tributary each (Big Beef, Chimacum, and Tahuya, respec- 
tively) where native populations had gone extinct. Two 
subpopulations, Dosewallips and Duckabush, were not 
deliberately supplemented. Supplementation programs 
were designed to reduce extinction risk and to speed recov- 
ery and recolonization while minimizing risks of deleteri- 
ous genetic, ecological, and demographic effects to 
supplemented and unsupplemented subpopulations 
(WDFW and PNPTT 2000; Tynan et al. 2003). Supplemen- 
tation hatcheries were scheduled to run for a maximum of 
three generations (12 years), and number of fry released 
into each tributary was based on available habitat. To best 
represent the demographics of each local donor population 
and to maintain existing genetic diversity, programs 
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followed protocols (Schroder and Ames 2004) in which 
supplementation broodstocks for each tributary were col- 
lected from local spawners randomly in proportion to the 
timing, weekly abundance, and duration of the total return 
from the same tributary where hatchery-propagated fish 
were to be released (each tributary had its own supplemen- 
tation program all with the goal of equal sex ratios and 
using no more than 50% of in-river spawners; see Table SI 
in Supplementary Information for broodstock details for 
each program). Matings employed partial factorial designs 
(Campton 2004; Busack and BCnudsen 2007) to maximize 
genotypic diversity and effective subpopulation size (Wa- 
ples and Do 1994; Withler and Beacham 1994). All hatch- 
ery fish were mass-marked either externally (fin-clipped, 
Quilcene program) or internally [all other programs 
received unique hatchery-specific otolith banding pattern 
marks (see Volk et al. 1987)]. Marks identified hatchery 
fish back to specific hatchery of origin from observed fin- 
clips or from otoliths collected from carcasses during 
annual spawner surveys (WDFW and PNPTT 2000; 
WDFW and PNPTC 2007). Because most hatchery fish 
were marked internally, after supplementation fish began 
to return, broodstock selected from returning spawners 
included unknown numbers of hatchery-born fish that 
were identified post-mortem. All supplementation program 
releases were mass-marked beginning with brood year 1997 
(and chum salmon generally return after 3 and 4 years), 
and from 2001 onward managers identified hatchery-born 
spawners from internal and external marks (Table 2). 
Although we have no direct measure of the reproductive 
success for hatchery-born spawners, the proportion of 
hatchery-born spawners was an indication of maximum 
expected hatchery influence. 

CoUections 

Collections consisted of summer chum salmon from HC 
and SJF subpopulations (Fig. 1; Table 1). Tissue samples 
for this study were collected from spawners in seven rivers 
and creeks in HC and two creeks in SJF from 1978 through 
2009 (Table 1). Samples collected prior to initiating sup- 
plementation programs were termed 'original' and samples 
collected after supplementation fish returned were divided 
by birth location, 'hatchery' or 'wild'. A few collections 
included samples from more than 1 year because there 
were too few samples per individual year (Table 1). The 
following summer chum salmon subpopulations had sup- 
plementation programs: Union, Lilliwaup, Hamma Ham- 
ma, Big Quilcene, Salmon, and Jimmycomelately (see 
Table 2 for numbers of hatchery-born fry planted per year 
and see Table S2 in Supplementary Information for years 
when hatchery-born spavmers returned). Although Big Beef 
Creek summer chum salmon were extirpated by 1990, we 



analyzed tissues (archived scales) collected in the creek 
prior to extirpation (1978, 1979) and after reintroduction 
with Quilcene stock (2004). Summer chum salmon in 
Dosewallips and Duckabush rivers (unsupplemented) were 
monitored for comparisons with supplemented subpopula- 
tions and to document and identify strays from supple- 
mentation programs. Spawners from all HC and SJF 
summer chum salmon subpopulations were censused 
annually for supplementation hatchery- and wild-born 
spawners in each river (post-mortem otolith analyses iden- 
tified hatchery-born spawners and their specific supple- 
mentation programs). Marked otoliths or fin-clip ratios 
from spawner carcasses estimated the number and propor- 
tion of hatchery-born spawners (Table 2), and scales were 
collected to age spawners as part of biological data for 
recovery monitoring (WDFW and PNPTC 2007). 

Genotyping 

We geno typed fish at 16 microsateUite loci [O/ce-3 (Buch- 
holz et al. 2001), Ofa'-l (Smith et al. 1998), Omy- 1011 
(Spies et al. 2005), One-lOl, One-102, One-106, One-108, 
One-Ill, One-lU (Olsen et al. 2000), One-18 (Scribner 
et al. 1996), Ots-l, Ots-2M, Ots-3,M (Banks et al. 1999), 
Ots-103 (Small et al. 1998), Ofs-G311 (Williamson et al. 
2002), and Ssfl-419 (Cairney et al. 2000)] for 2057 individ- 
uals from 43 collections (Table 1). DNA was extracted 
using silica membranes following the manufacturer's 
instructions (Macherey-Nagel). MicrosateUite loci were 
amplified in six multiplexes, and alleles were scored by two 
researchers prior to export for analyses (see Small et al. 
2009 for details of PCR, scoring, and binning). 

Statistical tests 

Genetic statistics were calculated to examine whether col- 
lections met expectations of random sampling and to iden- 
tify any problems with loci such as disequilibrium from 
large-allele dropout, null alleles, or scoring problems. 
Statistics assessed differences in genetic diversity between 
original samples, supplementation program samples 
(hatchery-born spawners and wild-born spawners collected 
in the same spawning areas), and samples from unsupple- 
mented tributaries and evaluated changes following supple- 
mentation using tests described below. The software 
COLONY2.0.0.1 (Wang 2008) calculated ftill-sibling family 
structure that might contribute to Ryman-Laikre effects 
through unequal representation of hatchery families among 
spawners. Samples were tested for departures from Hardy- 
Weinberg equilibrium (HWE) at each locus and across all 
loci using FSTAT 2.9.3 (Nei 1987; Goudet 2001) with 1000 
permutations. Departure from HWE can be an indication 
that samples contain family groups, or a strong year class, 
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Linkage (120 pairs) 



Region 


Tributary 


Year 


Name 


N 




Pvalue 


5% 


1% 


HC 


Big Beef 


1 978 


78BigB_0 


44 


—0.007 


0.6526 


8 


3 


HC 


Big Beef 


1 979 


79BigB_0 


43 


0.045 


0.0066 




3 


HC 


Big Beef 


2004 


04BigB_H 


45 


-0.003 


0.5768 


1 


4 


HC 


Union 


2000 


00Union_O 


54 


0.01 6 


0.1 570 




2 


HC 


Union 


2003 


03Union_W 


32 


-0.01 5 


0.7638 


5 


0 


HC 


Union 


2003 


03Union_H 


30 


0.013 


0.2706 


21 




HC 


Union 


2004 


04Union_H 


43 


0.014 


0.2255 


28 


m 


HC 


Union 


2008 


08Union_W 


47 


—0.022 


0.9165 




3 


HC 


Lilliwaup 


1985 


85Lilli_0 


23 


-0.037 


0.9450 




6 


HC 


Lilliwaup 


1 992 


92Lilli_0 


46 


0.007 


0.3291 




3 


HC 


Lilliwaup 


1 997-99 


97_99Lilli_W 


35 


—0.01 7 


0.8032 




H 


HC 


Lilliwaup 


2000-01 


00_01 LillLW 


35 


0.019 


0.1 504 




3 


HC 


Lilliwaup 


2001 


01 Lilli_H 


26 


—0.006 


0.591 5 




1 


HC 


Lilliwaup 


2002 


Q2Lilli_H 


47 


-0.020 


0.8725 




1 


HC 


Lilliwaup 


2003-04 


03_04Lilli_H 


53 


-0.006 


0.6558 


84 


1 


HC 


Lilliwaup 


2005-06 


Q5_Q6Lilli_W 


29 


-0.033 


0.9512 


26 


■ 


HC 


Lilliwaup 


2005 


Q5Lilli_H 


60 


0.000 


0.51 12 




2 


HC 


Lilliwaup 


2006 


Q6Lilli_H 


55 


0.022 


0.0698 


1 


2 


HC 


Lilliwaup 


2008 


Q8Lilli_H 


42 


0.007 


0.331 5 




3 


HC 


Lilliwaup 


2009 


Q9Lilli_H 


29 


0.025 


0.1052 


20 


4 


HC 


Hamma Hamma 


1999 


99Hamma_0 


34 


0.033 


0.0399 


22 


4 


HC 


Hamma Hamma 


2001 


01 Hamma_W 


53 


0.001 


0.4666 




2 


HC 


Hamma Hamma 


2003 


Q3Hamma_W 


40 


-0.013 


0.761 1 




2 


HC 


Hamma Hamma 


2001-03 


Q1_Q3Hamma_H 


27 


0.01 8 


0.1 830 




4 


HC 


Hamma Hamma 


2008 


Q8Hamma_W 


41 


0.000 


0.4835 




1 


HC 


Dosewallips 


1 992 


92Dose_0 


48 


—0.01 2 


0.7837 


6 


2 


HC 


Dosewallips 


2000 


*00Dose_W 


53 


0.029 


0.0369 


1 


2 


HC 


Dosewallips 


2003 


*03Dose_W 


46 


0.048 


0.0015 


6 


3 


HC 


Dosewallips 


2009 


*09Dose_W 


33 


-0.01 1 


0.7137 




0 


HC 


Duckabush 


1986 


86Duck_0 


57 


0.014 


0.1 713 




4 


HC 


Duckabush 


1992 


92Duck_0 


46 


-0.007 


0.6837 


9 


2 


HC 


Duckabush 


2000 


*00Duck_W 


46 


0.037 


0.0141 


13 


2 


HC 


Duckabush 


2003 


*03Duck_W 


47 


0.018 


0.1352 


9 


3 


HC 


Duckabush 


2009 


*09Duck_W 


34 


0.018 


0.1637 


1 


2 


HC 


Quilcene 


1 992 


yzi^uii_u 


50 


—0.01 2 


0.7824 


3 


1 


HC 


Quilcene 


2008-09 


no nn/^iiil \a/ 

Uo_uyQuii_w 


45 


0.009 


0.2346 


mm 

1 


4 


SJF 


Salmon 


1 986 


86Salmon_0 


42 


0.024 


0.6322 


1 


4 


SJF 


Salmon 


2000 


Q0Salmon_W 


25 


0.027 


u.oyjo 


1 


4 


SJF 


Salmon 


2000 


Q0Salmon_H 


32 


0.003 


0.1 133 


1 


3 


SJF 


Salmon 


2003-05 


Q3_Q5Salmon_W 


34 


0.005 


0.4324 


1 


■ 


SJF 


Salmon 


2003 


Q3Salmon_H 


27 


0.01 8 


0.4063 


■ 


■ 


SJF 


Salmon 


2004-05 


04 05Salmon H 


24 


0.131 


0.2010 


4 


0 


SJF 


Salmon 


2008 


08Salmon_W 


25 


-0.018 


0.0000 




2 


SJF 


Salmon 


2009 


09Salmon_W 


23 


0.014 


0.7968 




3 


SJF 


Jimmycomelately 


1986 


86JCL_0 


62 


-0.035 


0.2754 




■ 


SJF 


Jimmycomelately 


1 998-99 


98_99JCL_0 


16 


-0.019 


0.9910 




5 


SJF 


Jimmycomelately 


2000 


00JCL_O 


38 


0.019 


0.7256 






SJF 


Jimmycomelately 


2001 


01JCL_W 


58 


0.009 


0.1688 


45 


1 


SJF 


Jimmycomelately 


2001-04 


03_Q4JCL_H 


75 


-0.017 


0.2792 


80 




SJF 


Jimmycomelately 


2008-09 


08_Q9JCL_W 


33 


0.001 


0.8896 


11 


3 


SJF 


Jimmycomelately 


2008 


Q8JCL_H 


20 


-0.028 


0.4732 


12 


5 


SJF 


Jimmycomelately 


2009 


09JCL„H 


34 


-0.042 


0.8736 


1^ 


6 



Regions are abbreviated: Hood Canal = 'HC, Strait of Juan de Fuca = 'SJF'. Name abbreviations used throughout document are in 'Name' column. 
Name includes collection category (0 = original prior to supplementation program, H = Hatchery-origin, and W = Wild-origin during or after time of 
hatchery supplementation), and * indicates collections from tributaries that were not deliberately supplemented. Statistics include the Hardy-Wein- 
berg equilibrium value expressed by fis, and its associated Pvalue (underiined Fis values significant before Bonferroni corrections, and bold value was 
significant after Bonferroni corrections), the number of locus pairs (out of 1 20 pairs) in linkage disequilibria at the 5% and 1 % level (values in gray 
boxes were >B% of total pairs). Other genetic statistics, gene diversity and allelic richness, are presented in Table 4. 
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or more than one subpopulation, or related parents in the 
previous generation, or locus problems such as null alleles 
or selection. We tested whether genotypes at each locus 
were independent with the linkage disequilibrium permu- 
tation test in GENETIX 4.03 (Belkhir et al. 2001) using 500 
permutations. For each sample, GENETIX calculated the 
number of loci pairs in which 5% and 1% of the permuta- 
tions had a smaller value than the actual value for the sam- 
ple. Linkage disequilibrium can be an indication that the 
sample contains family groups or the offspring of matings 
of genetically distinct populations, or that the population is 
under selection, or that alleles have drifted due to a small 
subpopulation size. Basic genetic diversity measures includ- 
ing gene diversity (Nei's (1987) estimate of heterozygosity 
corrected to minimum eight individuals) and allelic rich- 
ness (average number of alleles per locus corrected to mini- 
mum eight individuals) were calculated using FSTAT. We 
used Wilcoxon signed-rank tests to test for significant dif- 
ferences among sample categories in gene diversity, allelic 
richness, and N,. that might signal effects associated with 
supplementation. Results for all tests were adjusted for 
multiple comparisons (sequential Bonferroni correction. 
Rice 1989) to an alpha level of 0.05. 

Two methods were used to estimate N^; N^. is a funda- 
mental parameter determining the evolutionary potential of 
a population, yet different techniques employ different 
assumptions and yield somewhat different results. Tracking 
changes in may indicate changes associated with supple- 
mentation (Antao et al. 2010; Hare et al. 2011). We esti- 
mated Ne using linkage disequilibrium (LD-Ne, Waples 
2006) in the program LDNe (Waples and Do 2008; setting 
the lowest frequency allele at 2% to avoid bias introduced 
by small collections) and estimated using a method 
based on maximum likelihood pairwise sibship analysis 
(SA-Ne, Wang 2009) implemented in the program COL- 
ONY 2.0.0.1. Because chum salmon have overlapping year 
classes, the values from these measures estimated the num- 
ber of breeders (Nj,) giving rise to the collection, rather than 
the actual N^. To calculate actual N^, the Nb is multiplied by 
the generation time, which is 3-5 years in chum salmon. 
Because generation time is variable in chum salmon (espe- 
cially variable in these recovering subpopulations), we sim- 
ply referred to the values calculated for a single-year mixed- 
age collection by LD and SA as and did not multiply by 
generation time. Because assumptions and results varied 
somewhat between methods, we calculated a harmonic 
mean (and variance) of values (HmNe) from the two 
methods, following Waples and Do (2010). We also 
calculated the ratio of HmNe to the census size for spawner 
types as a metric to investigate differences among sample 
types and changes over time (Hedrick 2005). 

Effective population size analyses were conducted from a 
metapopulation perspective (Gomez-Uchida et al. 2013): 
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we calculated the metapopulation before supplementa- 
tion for the earliest original program samples (00Union_O, 
85Lilli_0, 99Hamma_0, 92Quil_0, 86Salmon_0 and 
86JCL_0) and after supplementation for the latest collected 
wild-born samples (08Union_W, 05_06Lilli_W, 08Ham- 
ma_W, 08_09Quil_W, 09Salmon_W, and 08_09JCL_W). 
The metapopulation was first calculated according to 
the Stepping Stone model (Maruyama 1970) and the Island 
model (Wright 1943) using values from the LD-Ne and 
SA-Ng and then we calculated the harmonic mean of the 
metapopulation LD-N^ and SA-N^ for final metapopula- 
tion HmNe values (meta-HmNe) based on the Stepping 
Stone and Island models. For the Stepping Stone model the 
average immigration rate was estimated from the global 
_Fsx values for the two sets of six collections (original and 
wild-born), and the same global Fsr values were used in 
denominators for the Island model. 

Subpopulation comparisons 

A variety of measures assessed whether supplementation 
had induced changes in genetic attributes and genetic 
structure and whether strays from supplementation pro- 
grams into unsupplemented tributaries had altered popu- 
lation structure or substructure. To assess possible 
changes in genetic attributes (heterozygosity, allelic rich- 
ness, and HmNg), we used Wilcoxon sign rank tests to 
test for significant differences in these genetic attributes 
among collection categories; original, supplementation 
hatchery-born, and wild-born (see Table 1). We tested 
for significant differences in genotypic distributions 
among all samples using GENEPOP3.3 (Raymond and 
Rousset 1995), and examined temporal and spatial parti- 
tioning of pairwise genetic variance with pairwise F^r 
tests in GENETIX 4.03 (Belkhir et al. 2001). In pairwise 
-FsT tests we evaluated whether variance was significantly 
different from zero with 1000 permutations. A permuta- 
tion test implemented in FSTAT (comparison among 
groups of samples) was used to test for significant differ- 
ences in global fsx values among categories of samples 
(original, supplementation hatchery-born and wild-born, 
see Table 1), with 10 000 permutations. To assess 
changes associated with supplementation, we compared 
pairwise _Fst values among original samples and among 
wild-born samples collected after supplementation pro- 
grams had been initiated. These tests were conducted 
with and without samples from Dosewallips and Ducka- 
bush, the tributaries that were not deliberately supple- 
mented. To assess changes in variance associated with 
hatchery supplementation programs, we compared pair- 
wise _FsT values among samples of hatchery-born spaw- 
ners and wild-born spawners collected in the same 
tributaries. 



Genetic distances 

We assessed subpopulation structure using a principle 
coordinates analysis (PCoA) implemented in GenAlEx 6.5 
(Peakall and Smouse 2006; PeakaU and Smouse 2012). 
GenAlEx calculated pairwise _Fsx values among populations 
over all loci and conducted a PCoA of the pairwise values. 
The PCoA is a multivariate ordination technique that 
describes underlying patterns in a dataset. The PCoA gen- 
erates axes describing genetic variance in a dataset, here 
pairwise _Fsx values, with the first two axes usually describ- 
ing the maximum variance. As another means to view 
genetic relationships, we plotted CavaUi-Sforza and 
Edwards (1967) chord distances among collections in a 
dendrogram with 10 000 bootstrap replications using pro- 
grams within PHYLIP (Felsenstein 1993). 

Assignment tests 

To investigate whether fish sampled from the same tribu- 
tary were more or less likely to be genetically similar to 
each other following supplementation, we used self-assign- 
ment tests in GeneClass2 (Piry et al. 2004). GeneClass 
employs the Rannala and Mountain (1997) algorithm in a 
'leave one out' protocol and calculates the likelihood that 
an individual fish originated in the subpopulation in the 
tributary where it was sampled (home collection) based on 
the genotype of the fish and allele frequencies in the base- 
line collections, with the fish in question removed from its 
home collection in the baseline. There was no threshold 
likelihood value for assignment, the highest likelihood was 
accepted as the assignment and self-assignments were con- 
ducted on collections before and after supplementation. 
We tested for differences in percentage of assignments to 
home collection before and after supplementation using 
paired Student's f-tests. 

Isolation by distance 

Data were compared before and after supplementation for 
evidence of changes in isolation by distance (IBD) patterns 
(Slatkin 1993), using transformed pairwise f^x values [fsx/ 
(1--Fst)] as a genetic similarity measure (Rousset 1997). 
Geographical distances (kilometers) between mouths of 
streams were calculated using the most direct passage over 
open water. Mantel tests for association between pairwise 
_FsT values and distance, and reduced major axis regressions 
were performed using IBD 3.23 (Bohonak 2002). To assess 
changes that might have occurred as a result of supplemen- 
tation, we used ancova to compare IBD patterns among 
original collections to IBD patterns among collections 
of wild-born fish after supplementation began. We also 
compared IBD patterns of wild-born samples and 
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hatchery-born samples to evaluate differences between the 
sample groups. We conducted analyses with and without 
Dosewallips and Duckabush collections because strays from 
hatchery programs were documented in spawner surveys in 
those tributaries and may have supplemented those popu- 
lations. 

Results 

Supplementation and escapement 

Supplementation programs boosted abundance of spaw- 
ners and contributed from 5% to 96% of spawners to 
escapements in supplemented HC and SJF subpopula- 
tions (Table 2). The proportion of hatchery-born 
spawners decreased as successful programs terminated 
according to protocols (WDFW and PNPTT 2007). From 
2001-2011, strays from supplementation programs con- 
tributed from 0% to 33% (average = 10%) of spawners 
to Dosewallips and Duckabush (Table 2), which were not 
deliberately supplemented. Most strays to Duckabush and 
Dosewallips were from nearby supplementation programs 
in Hamma Hamma and Quilcene (WDFW and PNPTT 
2007). Similarly, most strays in other tributaries were 
from nearby supplementation programs and straying may 
be influenced by the natural exchange rate among these 
subpopulations. 

Subpopulation statistics 

Genotypic coverage averaged 95% and ranged from 84% 
(One-106) to 97% (Ofa'-l) per locus over all individuals 
{N = 2086). Because samples (tissues and archived scales) 
were from spawner carcasses of varying freshness, genotyp- 
ing success varied among samples, regardless of contempo- 
rary or historical status. In tests for Hardy-Weinberg 
equilibrium (HWE) at individual loci (Supplementary 
Information, Table S3), two tests out of 832 were signifi- 
cant after corrections for multiple tests and all samples 
except 08Salmon_W were in HWE in tests over all loci 
(Table 1). The Wilcoxon sign rank tests indicated no sig- 
nificant differences between original and wild-born samples 
in genetic diversity measures (gene diversity, allelic rich- 
ness) or HmNe (Fig. 2, Table 3, see Tables 1 and 4 for indi- 
vidual collection values). However, hatchery-born 
spawners had lower genetic diversity and significantly lower 
HmNg than wild-born spawners (differences in genetic 
diversity were not significant when Dosewallips and Ducka- 
bush were excluded from wild samples but differences in 
HmNe between wild- and hatchery-born samples remained 
significant). 

Linkage disequilibrium and pairwise sibship analyses 
suggested most collections included some family groups. 



ranging from a pair of siblings to large families (data not 
shown). There were 120 locus pairs examined per collec- 
tion, and most collections had six or more locus pairs in 
linkage disequilibrium at the 5% level and 12 collections 
had six or more locus pairs in linkage disequilibrium at the 
1% level (Table 1). The Wilcoxon signed-rank test indi- 
cated that hatchery-born spawners had significantly more 
locus pairs in linkage disequilibrium than wUd-born spaw- 
ners (P = 0.019). The collections from 02Lilli_H, 03_04Lil- 
li_H, and 03_04JCL_H were notable for high linkage 
disequilibrium: COLONY calculated a single family with 14 
fuU-siblings in 02Lilli_H (12% of pairwise relationships 
were fuU-siblings), seven families with between three and 
six full-siblings in 03_04Lilli_H, and six families with 
between three and eight full-siblings in 03_04JCL_H (4% 
of pairwise relationships were full-siblings). Siblings in col- 
lections were noted but not removed. 

Effective population size calculations 

We calculated and its 95% parametric confidence inter- 
val for each collection with two methods: linkage disequi- 
librium (LD-Ne) and pairwise sibship analysis (SA-Ne). In 
most calculations, the two values were similar and the 95% 
confidence intervals overlapped (see data for both methods 
for each collection in Supplementary Information Table S4 
'A11_N(,'). Differences arose for collections with high link- 
age and small samples sizes, which biased LD-ATe downward 
(Wang and Whitlock 2003; Waples and Gaggiotti 2006), 
and where data were mostly lacking for a single locus, 
which depressed SA-Ne in 79Big Beef sample. The har- 
monic means of values (HmNe) varied over time and 
space, both within and between tributaries (Table 4, 
Fig. 3). Where there were samples from hatchery-born and 
wild-born spawners collected in the same tributary in the 
same year, the hatchery-born samples generally had a smal- 
ler HmNe (e-g- 03Union_H and 03Union_W), except for 
samples from Lilliwaup Creek: in two comparisons of 
hatchery- and wild-born samples for that creek, the HmNe 
for the hatchery-born samples was equal to or greater than 
the HmNe for the wild-born sample. In the original Lilli- 
waup samples, the HmNe values were similar to those of 
Union River, but they declined to less than half by 1997. 
We lack wild-born samples from Lilliwaup Creek after 
2006, but in the most recent hatchery-born samples, the 
HmNe values were similar to original values, suggesting 
that diversity may be recovering in the Lilliwaup Creek 
subpopulation. The HmNe in the original Jimmycomelately 
Creek sample was roughly half the value of the most recent 
wild sample, suggesting that diversity is also recovering in 
the limmycomelately Creek subpopulation. The calculated 
ratios of HmNe to census size (HmNe/N in Table 4) were 
significantly lower in wild-born samples than in original 
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Figure 2 Plots of genetic statistics for samples grouped into original, hatchery, and wild categories. Genetic statistics are detailed in Tables 1 and 2. 



samples (Wilcoxon sign rank test, P= 0.0011 and 
P = 0.0137 with Dosewallips and Duckabush collections 
included and excluded, respectively), likely reflecting 
increases in census size (N) throughout the restoration pro- 
gram. Further, in an anova the ratio of HmNe to census 
size for the wild-born fish was negatively correlated with 
the number of years of supplementation (fi^es = 11-2, 
P = 0.001). 

The global Fst values for the six original and six wQd- 
born samples were 0.029 and 0.021, respectively 
(P = 0.2414), resulting in calculated immigration rates of 
8.37 and 11.65, respectively. The meta-HmNg for the origi- 
nal and wild-born samples was 488 and 591, respectively, 
with the Stepping Stone model, and 502 and 603, respec- 
tively, with the Island model (see Supplemental Informa- 



tion Table S4 for LD-Ne and SA-Ne values of six original 
and six wild-born samples and calculations for meta-HmZVe 
values). Paired f- tests for LD-Ne and SA-Ne values indi- 
cated no significant differences between values for original 
and wild-born fish (P = 0.1873 for LD-Ne values for origi- 
nal versus wild-born fish, P = 0.3697 for SA-N^ values for 
original versus wild-born fish). Although meta-HmZVe 
increased following supplementation, the difference was 
not significant. 

In one further consideration of N^, we computed the 
harmonic mean of the calculated values over all collec- 
tions with hatchery-born and wild-born spawners com- 
bined in collections and separated (see Supplementary 
Information, Table S5 'All_HmNe'). The calculated value 
for the wild-born spawners (78.95 ± 33.86) was larger 
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Table 3. Table of averages (avg) and harmonic means (hmean) of genetic statistics for categories of samples (original = O, hatchery-born = H, wild- 
born = W) detailed in Tables 1 and 4, and P values for Wilcoxon signed-rank tests for comparisons 
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Analyses were conducted with and without Duckabush and Dosewallips samples (DoseDuck). Gene diversity is Nei's estimate of heterozygosity cor- 
rected to a sample size of 6, and allelic richness is average number of alleles per locus corrected to a samples size of 6. 



than the value calculated for the hatchery- and wild-born 
spawners combined (55.19 ± 37.44), but the confidence 
intervals overlapped. Also, within some single-year examin- 
ations, the ratios of HmNe/N were higher in uncombined 
than in combined collections, suggesting Ryman-Laikre 
effects from combining hatchery- and wild-born fish into 
single collections. 

Genetic variance patterns within and between 
subpopulations 

Pairwise genotypic and fsx tests were mostly congruent 
(see Supplementary Information Table S6 for pairwise test 
values and their associated P values). Pairwise tests indi- 
cated temporal stability within most subpopulations, with 
the exception of LiUiwaup Creek (see Fig. 3 for a plot of 
pairwise Fst values along principle coordinate axes). 
Among LiUiwaup Creek samples, 37 of 66 pairwise geno- 
typic comparisons were significant and 34 of 66 pairwise 
comparisons were significant. Union River samples also 
differed from all other HC samples and are the most iso- 
lated geographically (Fig. 1). There was little differentiation 
among samples from Dosewallips, Duckabush, Hamma 
Hamma, and Quilcene rivers. The SJF samples differed 
from each other and were distinct from HC samples. The 
original 1970s Big Beef Creek samples were most similar to 
those of Union River and LiUiwaup Creek, and the reintro- 
duced sample was most similar to the sample from Quil- 
cene River, its broodstock source, and consequently similar 
to DosewaUips, Duckabush, and Hamma Hamma rivers 
samples (Fig. 3). 

Global Fsi values for different categories of samples did 
not change significantly throughout the program (Table 5). 
Although the global Fgr value for wild-born samples was 
lower than the value for original samples, the decrease fol- 



lowing supplementation was not significant. The global Fsr 
for original samples matched the global Fsr for the hatch- 
ery samples, suggesting that the hatchery programs cap- 
tured similar genetic variance among subpopulations. 
Including original and wild-born coUections from Dose- 
waUips and Duckabush in computations lowered global f si- 
values, but changes were still not significant. 

Genetic clusters identified in PCoA and dendrogram 

Population centers formed two major clusters in the princi- 
ple coordinates analysis plot (Fig. 3), and there was no dis- 
tinction detected between hatchery and wUd samples from 
the same tributaries. The first axis explained 45% of the 
genetic variance: SJF samples clustered on the left, and HC 
samples clustered on the right. The second axis explained 
12% of the variance: The HC populations divided loosely 
into two to three clusters along the second axis. The lowest 
cluster included samples from DosewaUips, Duckabush, 
Hamma Hamma, and Quilcene rivers as well as the 2004 
coUection from Big Beef Creek. The topmost cluster 
included the original samples from Big Beef Creek (these 
separated along the third axis; the 1979 coUection was miss- 
ing most of one locus and plotted distantly from the 1978 
coUection), Union River samples, and the older samples 
from LiUiwaup Creek. The more recent LiUiwaup samples 
clustered mainly in the center, supporting high genetic drift 
in this population. The placement of the original Big Beef 
samples indicated that before reintroduction, summer 
chum salmon in Big Beef Creek were genetically more simi- 
lar to other subpopulations from low-elevation tributaries 
on the east side of HC rather than the geographically closer 
subpopulations in tributaries originating in the Olympic 
Mountains (Dosewallips and Duckabush) on the west side 
of HC. A dendrogram of genetic distances among samples 
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Table 4. Population statistics for summer chum salmon collections including gene diversity ('GeneDiv', Nei's estimate of heterozygosity corrected to 
a sample size of 8) and allelic richness {'Ar', average number of alleles per locus corrected to a sample size of 8) 
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01_03Hamma_H 


0.8151 


7.42 


71 


14 


390 


0.18 


HC 


08Hamma_W 


0.8127 


7.56 


163 


33 


1371 


0.12 


HC 


92Dose_0 


0.8177 


7.54 


151 


22 


655 


0.23 


HC 


OODose_W 


0.8162 


7.50 


91 


11 


1260 


0.07 


HC 


03Dose_W 


0.8178 


7.53 


120 


24 


6510 


0.02 


HC 


09Dose_W 


0.8167 


7.59 


106 


19 


1093 


0.10 


HC 


86Duck_0 


0.8199 


7.46 


101 


10 


234 


0.43 


HC 


92Duck_0 


0.8169 


7.46 


155 


22 


617 


0.25 


HC 


OODuck_W 


0.8121 


7.32 


100 


13 


464 


0.22 


HC 


03Duck_W 


0.8026 


7.27 


92 


15 


1600 


0.06 


HC 


09Duck_W 


0.8062 


7.34 


133 


24 


2499 


0.05 


HC 


92Quil_0 


0.8063 


7.35 


125 


16 


743 


0.17 


HC 


08_09Quil_W 


0.8109 


7.24 


71 


9 


5353 


0.01 


SJF 


86Salmon_0 


0.7842 


6.67 


43 


6 


582 


0.07 


SJF 


OOSalmon_W 


0.8116 


7.06 


133 


64 


439 


0.30 


SJF 


OOSalmon_H 


0.8069 


6.74 


64 


13 


407 


0.16 


SJF 


03_05Salmon_W 


0.8173 


7.31 


90 


13 


7642 


0.01 


SJF 


03Salmon_H 


0.7965 


6.70 


36 


9 


1866 


0.02 


SJF 


04_05Salmon_H 


0.8134 


6.77 


67 


12 


4203 


0.02 


SJF 


08Salmon_W 


0.8239 


7.33 


111 


25 


1544 


0.07 


SJF 


09Salmon_W 


0.7991 


6.92 


88 


28 


1218 


0.07 


SJF 


86JCL_0 


0.7776 


6.23 


32 


4 


292 


0.11 


SJF 


98_99JCL_0 


0.7593 


5.99 


21 


7 


105 


0.20 


SJF 


00JCL_O 


0.7589 


6.00 


22 


6 


55 


0.40 


SJF 


01JCL_W 


0.7795 


6.27 


36 


4 


251 


0.14 


SJF 


03_04JCL_H 


0.7665 


6.02 


25 


5 


1427 


0.02 


SJF 


08_09JCL_W 


0.7831 


6.46 


65 


11 


1099 


0.06 


SJF 


08JCL_H 


0.7936 


6.67 


56 


15 


481 


0.12 


SJF 


09JCL_H 


0.7638 


6.21 


46 


10 


2102 


0.02 



Harmonic mean Ne (HmWj.) was calculated for We's from linkage disequilibrium (LD-We) and pairwise sibship analysis (SA-N^). Escapement (Escape) 
was calculated using area under the curve and fish counts at traps and sums wild- and supplementation hatchery-origin escapements to natural 
spawning grounds. The HmWe/N is the ratio of HmW^ to escapement. 
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Figure 3 Principle coordinates plot of pairwise Fst values among summer chum salmon collections from Hood Canal and Strait of Juan de Fuca. 
Name abbreviations follow Table 1 and colors for categories follow Fig. 2, with the exception of Big Beef (in black): '78 and 79 Big Beef collections 
were sampled prior to extinction, and the '04 Big Beef collection was derived from Quilcene (Quil) broodstock. 



(Fig. 4) displayed similar overall structure. The SJF and 
HC clusters separated with 100% bootstrap support, with 
JCL forming a supported sub-branch. The Dosewallips, 
Duckabush, Hamma Hamma, and Quilcene samples 
formed an unsupported cluster, and the original samples 
from Lilliwaup Creek clustered with 79% bootstrap sup- 
port on a branch with Union River and the original Big 
Beef Creek samples. The remaining Lilliwaup Creek sam- 
ples clustered on three unsupported branches. There were 
no distinctions evident between hatchery-born and wild- 
born samples collected in the same tributaries. 



to a closely related population. In presupplementation col- 
lections, self-assignments were low to indistinguishable 
from random in Dosewallips, Duckabush, Hamma Ham- 
ma, and Quilcene rivers. If we combined results from these 
rivers, assignments back to the combined 'HoodMet' were 
high (Table 6). Self-assignments were relatively high for 
Salmon and Jimmycomelately creeks and moderate in 
Union, Lilliwaup, and Big Beef creeks. Following supple- 
mentation, self-assignments decreased in all but Union 
River and Hamma Hamma River, but the decrease was not 
significant (Student's paired f-test, P = 0.054). 



Assignment test 

Assignment tests supported a metapopulation structure 
with straying among close subpopulations, especially along 
the west side of HC (Dosewallips, Duckabush, Hamma 
Hamma, and Quilcene, Table 6). Where high gene flow 
occurs, an individual may be equally likely to have been 
sampled in two or more baseline populations or may assign 



Isolation by distance 

We conducted IBD analyses for all categories of samples 
(Fig. 5A-D), both including and excluding Dosewallips 
and Duckabush samples. In all analyses, physical distance 
explained a significant amount (all > 70%) of the genetic 
variance among samples (all P < 0.001). In comparisons 
among sample categories, because the tests for homogeneity 
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Table 5. Table of global F^j values for categories of samples and P val- 
ues for permutation test comparison among category values 





W samples 


Global fsT 


Pvalue 


Original, all samples 


12 


0.024 


0.294 


Wild, all samples 


21 


0.020 




Original, no Dose no Duck 


9 


0.032 


0.122 


Wild, no Dose no Duck 


15 


0.021 




Original, only Dose and Duck 


3 


0.001 


0.522 


Wild, only Dose and Duck 


6 


0.003 




Hatchery 


16 


0.032 


0.098 


Wild, no Dose no Duck 


15 


0.021 





The 'N samples' are the number of samples in the category. Big Beef 
Creek samples were excluded because the program was a reintroduc- 
tion. 

among regressions in the ancova indicated significant dif- 
ferences between regressions, we were unable to compare 
regressions statistically. However, the IBD slope was steeper 
for comparisons among original samples (Fig. 5A,B) than 
the slope for comparisons among wild samples, suggesting 
that genetic variance among populations decreased with 
supplementation. The IBD slopes were similar for compari- 
sons among hatchery-born and wild-born samples 
(Fig. 5C), but genetic distances were greater among hatch- 
ery-born fish. When analyses were limited to comparisons 
involving Dosewallips and Duckabush samples (Fig. 5D), 
the slope for the original samples was slightly steeper than 
the slope for the wild samples. 

Discussion 

This study joins a growing body of research on the genetic 
effects of supplementation hatchery programs on wild fish 
populations. We examined a time series of a summer chum 
salmon metapopulation (a designated ESU) from before, 
during, and after supplementation, providing a temporal 
perspective of supplementation and responses by threa- 
tened subpopulations. We found that after several years of 
hatchery supplementation, there was little change in genetic 
diversity and harmonic means of effective population sizes 
(HmNe) in wild-born salmon returning to spawning areas 
throughout the ESU. However, genetic distances within the 
metapopulation decreased and assignments back to collec- 
tion of origin decreased following supplementation, sug- 
gesting higher gene flow within the metapopulation and 
lower genetic drift in subpopulations. Hatchery-born spaw- 
ners collected in the same locations as wild-born spawners 
usually had lower genetic diversity and smaller HmNe val- 
ues, suggesting that supplementation hatchery programs 
sampled a subset of genetic diversity in the target subpopu- 
lations and that factorial matings of hatchery broodstocks 
may have decreased HmNg. Higher linkage disequilibrium 
in hatchery-born fish indicated potential for Ryman-Laikre 




99Hamma_W 
03 OAUWU» 



Figure 4 Neighbor-joining dendrogram of Cavalli-Sforza and Edwards 
chord distances among samples. Bootstrap values over 70% are plotted 
on the tree nodes. 



effect, and combining hatchery- and wild-born fish into 
single collections lowered ratios of effective population 
sizes to census sizes. However, with the exception of the 
Lilliwaup Creek subpopulation, suspected Ryman-Laikre 
effects did not increase genetic drift in wild-born spawners. 
Thus, possible negative impacts of supplementation 
appeared minimal and likely diminished as wild-born 
spawner abundance increased, suggesting that population 
recovery in this summer chum salmon ESU was unim- 
paired by supplementation. 

Supplementation and conservation 

In the Pacific Northwest, supplementation hatcheries using 
local broodstocks have received increased interest as a tool 
in fish conservation to preserve native genetic diversity 
and maintain a foundation of response to environmental 
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Figure 5 Plots of isolation by distance patterns (transformed pairwise F^j value (fsiAl - Fst) versus geographical distance) for different categories of 
samples: plot (A) all original samples and all wild samples, plot (B) original and wild samples (minus Dosewallips and Duckabush), plot (C) hatchery 
and wild samples (minus Dosewallips and Duckabush), and plot (D) only original and wild samples in comparison with Dosewallips and Duckabush. 
Plot (A) includes data from plots (B and D). 



variation. Traditional fish hatcheries were a mainstay of 
fish management and harvest augmentation, but they 
introduced problems when nonlocal, hatchery-adapted fish 
interacted with wild fish (Araki et al. 2008; Christie et al. 
2012). Supplementation hatcheries were designed to reduce 
problems associated with traditional hatcheries in that 
domestication selection was minimized by using in-river 
broodstocks originally composed of wild-origin fish (subse- 
quent broodstocks would be mixtures of hatchery- and 
wild-origin fish) and by reforming hatchery practices. 
Hatchery fish from the same gene pool as wild fish were 
expected to be similarly adapted to the local environment. 
Here we discuss studies documenting impacts from supple- 
mentation programs and hypothesizing mechanisms lead- 
ing to differences between hatchery and wild fish in 
relation to chum salmon. 

One selective force in supplementation hatchery pro- 
grams is juvenile residence in a hatchery environment 
(McClure et al. 2008). Hatchery effects may be relatively 
lower for chum salmon because they out-migrate shortly 
after emergence, spending minimal time under hatchery 
rearing conditions. In contrast, most supplementation 
studies have been conducted in species with extended juve- 
nile freshwater residence times [steelhead, Atlantic salmon 
{Salmo salar), coho and Chinook salmon (reviewed by 
Araki et al. 2008 and Chilcote et al. 2011)]. In these 



species, juvenile residence time in freshwater is a critical 
component of their life history, and unless they are released 
as unfed fry, hatchery-origin juveniles often spend a mini- 
mum of 1 year in a hatchery environment before out- 
migrating. Kostow (2004) observed differences in steelhead 
juvenile phenotypes and lower survival in hatchery-origin 
juveniles from the same parent pool. In a common garden 
experiment, Chittenden et al. (2011) found that for coho 
salmon juvenile rearing environment was a key influence 
on smolt size and multiple behaviors, regardless of ances- 
try. Theriault et al. (2011) implicated sexual selection at 
the juvenile stage, especially on males, as a factor decreasing 
relative reproductive success (RRS) in hatchery-origin coho 
salmon spawning in the wild. 

Another selective force in supplementation hatchery pro- 
grams is human-controUed spawning (McClure et al. 
2008), and chum salmon would be impacted similar to 
other species. Artificial spawning bypasses mate selection, 
and Consuegra and Garcia de Leaniz (2008) found a posi- 
tive association between disassortative mating (natural 
spawners choose mates with dissimilar MHC profiles) and 
parasite resistance: In Atlantic salmon offspring parasite 
load decreased their RRS. Further, the summer chum sal- 
mon supplementation programs used factorial matings, 
which likely decreased effective population size relative to 
single-pair matings (Abadia-Cardoso et al. 2013). Multiple 
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factors likely contribute to the phenomenon that hatchery- 
origin fish and wild-origin fish with hatchery ancestry are 
less productive under natural spawning conditions (Araki 
et al. 2008; Chilcote et al. 2011). 

In other studies, supplementation impacts were unde- 
tected or equivocal. In Berejikian et al.'s (2009) study, 
hatchery- and wild-born summer chum salmon spawners 
(Quilcene fish that were included in this study) mated ran- 
domly and there was no significant difi«rence in ofi^pring 
produced per spawner type. Sharpe et al. (2010) found no 
significant difference in reproductive success of wild-origin 
and first-generation hatchery-origin steelhead derived from 
the same gene pool. Further, genetic diversity in steelhead 
populations supplemented with native broodstock 
remained unchanged after 20 years (Heggenes et al. 2006) 
and 58 years (Gow et al. 2011). There was also no loss in 
genetic diversity (Eldridge and Killebrew 2008) or fitness 
(Hess et al. 2012) in Chinook salmon supplemented with 
native broodstock. In this study, we lack parentage data to 
address RRS (see Christie et al. 2012) and address only 
changes in genetic parameters. Yet, coupled with the study 
by Berejikian et al. (2009), the increases in wild-born spaw- 
ner abundance and increases in recruits per spawner sug- 
gested that hatchery-born HC summer chum salmon 
spawned successfully, which may have contributed to 
increased census sizes and lower extinction risks and con- 
tributed to recovery in the metapopulation (Wang and 
Ryman 2001). The lower HmNJN ratio in wild-born spaw- 
ners could also be an indication that some of their parents 
were hatchery-born fish that, although abundant, had 
lower genetic diversity than wild-born fish because of 
unequal hatchery family sizes. 

Metapopulation structure and supplementation 

Contemporary genetic structure in HC summer chum sal- 
mon follows an IBD pattern, similar to original genetic pat- 
terns (Phelps et al. 1994; Small et al. 2009), in which the 
amount of genetic exchange depended on distance between 
spawner groups (Schtickzelle and Quinn 2007). Interest- 
ingly, the original collections from Big Beef Creek were 
most similar genetically to the collection from Union River, 
which is roughly 80 km distant, rather than to the closest 
tributary, Dosewallips River, roughly 8 km west across HC 
(see Fig. 1). This supports an ecoregional association 
between summer chum salmon inhabiting lowland streams 
on the Kitsap Peninsula on the east side of HC as suggested 
by Sands et al. (2009). 

The relationship between genetic distance and geograph- 
ical distance changed following supplementation in this 
summer chum salmon metapopulation - gene flow 
increased and genetic distances decreased - but the basic 
IBD pattern remained consistent over time. In another sup- 
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plementation program for coho salmon, genetic distances 
increased following supplementation, which was attributed 
to genetic drift, bottlenecks, and varying success of 
between-river transfers (Eldridge et al. 2009). In HC sum- 
mer chum salmon, each supplementation program used 
river-specific broodstocks rather than a common brood- 
stock for all rivers. Thus, between-river transfers would 
occur via strays (nonhoming), which are documented in 
WDFW and PNPTC (2007). As a byproduct of increasing 
abundance, supplementation may have increased straying 
and decreased differentiation or perhaps restored straying 
to levels that existed prior to supplementation when sub- 
populations were more abundant. The pairwise _Fst values 
and assignment tests support that subpopulations are well 
connected through gene flow, which may have contributed 
to the increases in HmNe values. Alternatively, prior to 
supplementation census sizes were low and genetic drift 
may have increased differentiation; following supplementa- 
tion and the fishing moratorium, increased census sizes 
may have countered genetic drift and decreased differentia- 
tion. 

In summary, our results show that, although changes 
were not significant, in the HC summer chum salmon ESU 
contemporary genetic diversity and effective population 
sizes increased from original levels, suggesting that as natu- 
ral production increased, possible negative impacts from 
supplementation (such as depressed from Ryman-Lai- 
kre effects) diminished. We suspect that negative genetic 
impacts may be less or shorter-lived for chum salmon 
because their juveniles spend little time in the hatchery 
environment. We wiU continue to monitor this metapopu- 
lation to assess the long-term impacts of supplementation 
and success of recovery efforts. 
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